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Abstract 

We show that the Sudakov factor from the resummation of double loga- 
rithms ln(s/fcy) contained in the distribution functions is responsible for the 
&t smearing mechanism employed in the next-to-leading-order QCD («a s 2 ) 
calculations of direct photon production, s is the center-of-mass energy, and 
kr the transverse momentum carried by a parton in a colliding hadron. This 
factor exhibits the appropriate s-dependent Gaussian width in such that 
our predictions are in good agreement with experimental data. 
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I. INTRODUCTION 



It was pointed out some time ago QXJ] that the global QCD analysis of 
the direct photon production processes from both fixed-target and collider 
experiments @, [|, [|, |j, || [7], || || has led to a puzzle: Most data sets show 
a steeper pt distribution than the next-to-leading-order (NLO) QCD (acr 2 ) 
predictions, pt being the transverse momentum of the direct photon. This 
behavior can not be explained by either global fits with new parton distribu- 
tion functions or improved photon fragmentation functions. To resolve the 
puzzle, a Gaussian type broadening of the transverse momentum carried 
by initial-state partons in colliding hadrons has been introduced [1]]. This 
smearing function enhances the low end region of the spectrum in p^ more 
than the high end, such that the QCD predictions have a steeper distribu- 
tion. The recent publication of the E706 direct photon data set [|l(| has 
confirmed the evidence for the kx effect. However, the physical origin of this 
effect and of its center-of-mass-energy dependent Gaussian width has not 
been understood yet. 

In this paper we shall propose the mechanism that is responsible for the 
A/p-smearing effect in direct photon production. It has been observed that 
large double logarithms ln 2 (s/A;|,) are generated, s being the center-of-mass 
energy, when the transverse degrees of freedom of the partons are taken into 



account ||11|| . The large logarithms are absorbed into parton distribution 



functions, and their all-order summation leads to a Sudakov factor, which 
describes the perturbative distribution of the partons in kx for different s. 
Since the Sudakov factor gives strong suppression at large k T , it resembles a 
Gaussian function with s-dependent width. These characteristics are qual- 
itatively consistent with those of the smearing effect mentioned above. We 
shall demonstrate quantitatively that the inclusion of the Sudakov factor in- 
deed modifies the NLO QCD («a s 2 ) predictions in the desired way, and the 
data of direct photon production can be explained. 

Due to the transverse degrees of freedom of the partons, the resummation 
of double logarithms should be performed in the impact parameter b space 
[ p| , |T2f , which is conjugate to kx- Hence, it is ln 2 (s6 2 ) that are resummed. 
After deriving the Sudakov factor, we Fourier transform it back to the kx 
space, and convolute it with the naive NLO factorization formula for direct 
photon production. Note that an analytical expression of the transformed 
Sudakov factor does not exist, and must be obtained numerically. However, 
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to make its smearing effect transparent, the Sudakov factor is parametrized 
by a Gaussian function exp(— T 2 b 2 /4) with the width T, which corresponds 
to exp(— fcy/r 2 ) in the k T space. We find that the s dependence of T is 
consistent with that employed in |l[], implying the success of our analysis. 

In Sect. II we apply the Collins-Soper-Sterman resummation technique 
P~2]j to direct photon production from hadron collisions, and present the 



explicit expression of the Sudakov factor. In Sect. Ill the Sudakov factor 
is parametrized as a Gaussian function, and the s dependence of its width 
is examined. The numerical results are compared with the data of direct 
photon production. Section IV is the conclusion. 



II. FACTORIZATION AND RESUMMATION 

In this section we derive the factorization formula for direct photon pro- 
duction 

h(pi) + h(p 2 )^ 1 (p T )+X , (1) 

in the collision of the hadrons h, and resum the involved large double log- 
arithms into a Sudakov factor. The hadron momenta are assigned as p± = 

(Pi,0,Ot) and p 2 = (0,^2,0^) with pf = P2 = px is the transverse 

momentum of the direct photon that is measured. The tree-level diagrams 
are shown in Figs. 1(a) and 1(b), in which both quarks and gluons out of 
the hadrons contribute. The partons carry the momenta £iPi + kjr, where 
i — 1, 2, are the longitudinal momentum fractions and kjr are the transverse 
momenta, kjr flow into the direct photon, and make up the momentum p?. 

We then consider higher-order corrections to Fig. 1(a) from Figs. 1(c)- 
1(h). The discussion for Fig. 1(b) is similar. Fig. 1(c), the self-energy cor- 
rection to a parton, and Fig. 1(d), the loop correction with a real gluon 
connecting the two partons from the same hadron, contain both collinear di- 
vergences from the loop momentum / parallel to p, and soft divergences from 
small I. Since the soft divergences cancel asymptotically as shown in ||11|| , the 
double logarithmic corrections are mainly collinear. Therefore, Figs. 1(c) and 
1(d) are absorbed into a distribution function associated with the hadron 
i. In most of kinematic regions the self-energy correction to the outgoing jet 
in Fig. 1(e) does not give collinear divergences. Because of the cancellation 
of soft divergences, Fig. 1(e) is absorbed into a hard scattering amplitude 
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H, which corresponds to the parton- level differential cross section da/dpx- 
Similarly, other radiative corrections to the box diagram are also infrared 
finite and grouped into H . 

The absorption of the irreducible diagram in Fig. 1(f) is more delicate. It 
contains collinear divergences from I parallel to p^ and soft divergences. For I 
parallel to p±, we replace the parton line in hadron 2 by an eikonal line in the 
direction n — v% — p^j ' J~s/2 on the light cone, and factorize the gluon into 
the distribution function of hadron 1 from the full scattering amplitude. The 
replacement by an eikonal line holds for both quarks and gluons as verified 
in The treatment of the gluon with I parallel to P2 is the same, but 

the direction of the eikonal line is n — V\ — p\j \J s/2. The factorization of 
the vertex-correction diagrams in Figs. 1(g) and 1(h) is similar to that of 
Fig. 1(f). In the leading region with the loop momentum I parallel to pi, 
they are assigned to the corresponding distribution functions. Certainly, soft 
divergences in the above diagrams cancel asymptotically. At last, when I 
is hard, Figs. l(c)-l(h) are grouped into the hard scattering amplitude H 
defined above. 

The eikonal lines on the light cone introduced above, collecting collinear 
gluons, are essential for the factorization of irreducible diagrams in Figs. 1(e)- 
1(h), and for the gauge invariance of the parton distribution functions. How- 
ever, to implement the resummation technique, we allow the gauge vector 
n to vary away from the light cone (n 2 ^ 0) temporarily. It will be shown 
that the Sudakov factor turns out to be n- independent. After completing 
the resummation, n is brought back to the light cone, and the gauge invari- 
ance of the parton distribution functions is restored. That is, the arbitrary n 
appears only at the intermediate stage of the formalism and as an auxiliary 
tool of the resummation. 

When the transverse degrees of freedom of the partons are taken into ac- 
count, the factorization must be performed in the b space, which is conjugate 
to k T ]nj. Hence, we arrive at the general factorization picture for direct 



photon production shown in Fig. 2, and the corresponding formula, 

x%,6. s, b, n) exp(zp T ■ b) , (2) 
/j, being a renormalization and factorization scale. We have assumed the sin- 
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gle b dependence in the above expression in order to compare our formalism 
with the naive kx smearing employed in the literature. A rigorous factor- 
ization formula for direct photon production from hadron collisions will be 
presented elsewhere. 

In the axial gauge n- A = the eikonal line disappears, and only Figs. 1(c) 
and 1(d) give double logarithmic corrections to the distribution functions 
ll|]. We demonstrate how to resum the double logarithms from Figs. 1(c) 
and 1(d) by considering only the hadron 1. The resummation for the hadron 
2 is the same. The essential step in the resummation is to derive a differential 



equation pid(f)/dpx = C(p [O]], where the coefficient function C contains only 
single logarithms, and can be treated by RG methods. In the axial gauge n 
goes into the gluon propagator, (—i/l 2 )N flu (l), with 

N , (l)=r ±_ + „._"_. (3) 

Because of the scale invariance of N^ v in n, <p must depend on pi (or s) 
through the ratio (pi ■ n) 2 /n 2 , implying that the differential operator djdpX 
can be replaced by d/dn using a chain rule, 

d ~ n 2 d ~ 

PiTT^ = v la —<f) . (4) 

dpj Vi ■ n dn a 

The operator d/dn a applying to N^ v gives 

n 2 d 

via-j—N"" = v a {N^ a l v + N au F) , (5) 



v% ■ n dn, 
with the special vertex 



n 2 vi a 



, (6) 

The momentum Z M (l v ) appearing at the end of the differentiated gluon line 
is contracted with a vertex the gluon attaches, which is then replaced by 
the special vertex. The contraction of hints the application of the Ward 
identity. Summing all the diagrams with different differentiated gluons, the 
special vertex moves to the outer end of the parton line due to the Ward 
identity. We obtain the derivative, 

^-40 = 20', (7) 
dpj 
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described by Fig. 3(a), where the square in the new function 0' represents 
the special vertex v a . The coefficient 2 comes from the equality of the two 
new functions with the special vertex on either of the two parton lines. 

The collinear region of the loop momentum I is now not important be- 
cause of the factor \/{n-l) in v a with nonvanishing n 2 . Therefore, the leading 
regions of I are soft and hard, in which the subdiagram containing the spe- 
cial vertex can be factorized from 0' into a function K and a function G, 
respectively. The remaining part is the original distribution function 0. The 
differential equation is then expressed as 

, d ~ 



Pi 



dpi 



= 2 K{bn,a 8 {n)) + G(pf/fi,a s (fi)) 0, 



(8) 



with the general diagrams of K and G shown in Fig. 3(b). The sum K + G 
is exactly the coefficient function C mentioned above. It has been made 
explicit that K contains the single small scale 1/6 and G contains the single 
large scale pf. 

The 0(a s ) contributions to K from Fig. 3(c) and to G from Fig. 3(d) are 
given by 

1 

( i*l 



K 



G 



ig Cji 1 



d 4 ~ e l 
d 4 ' 



4-e 



£ + 2vr^(/ 2 )e ilT ' b 



(27T) 




5K 



5G 



(9) 
(10) 



Xp + i) 2 

5K and 5G being the corresponding additive counterterms. The color factor 
C is Cp{= 4/3) if the parton is a quark, and N c (= 3) if the parton is a gluon 
[ |l"3|j . The factor e llr ' h in Eq. (|9|), which is associated with the second diagram 
(real gluon emission) in Fig. 3(d) JTT], [12|], renders K free of infrared poles. 
That is, 1/6 serves as an infrared cutoff of the loop integral. Note that K 
vanishes in the asymptotic region 6 — > 0, reflecting the soft cancellation as 
stated before. The second term in Eq. (|10|), as a soft subtraction, ensures 
that the momentum flow in G is hard. 

A straightforward calculation of Eqs. (0) and (|10D gives the pole terms 
5K = —5G. For details of the calculation, refer to [[l l|]. Since K and G 
contain only single soft and ultraviolet logarithms, respectively, they are 
treated by RG methods: 

i>^r(: ■ (11) 



d K 
d\i 



d[i 
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The anomalous dimension of K, Xk = fidSK/d/j,, is given, up to two loops, 



by m 



A 



K 



C + 



71 



c 




7"/ 



(12) 



with rif the number of quark flavors, and Ca = 3 a color factor. In solving 
Eq. (|TTD , we make the scale /i evolve to the infrared cutoff 1/6 in K and to 
p + in G. The RG solution of K + G is written as 

" dp, 



K(bfi, a s (fi)) + G(p + /fi, a 8 (p)) 

H/b n 

Substituting Eq. ( |T3"D into (§), we obtain the solution 



4>(€i,Pi,b,n) = f(€ipt,b)4>(£i,b,fj,) 
with the Sudakov factor 

f^pf dp 

/(CiPi , 6) = exp ° 



(13) 



(14) 



(15) 

'1/6 p Jl/fe /X 

We have set the lower bound of the variable p to 1/6, and the upper bound 
to £iPi". Note that Eq. (|15|) is defined only for ^pf > 1/6. For £iPi < 1/6, 
we require / to be equal to unity. We also set / to 1 as f > 1 (T^j, which 
occurs in the small 6 region. The radiative corrections in this short-distance 
region should be absorbed into the hard scattering amplitude H, instead of 
into the distribution function, giving its Sudakov evolution. The wave func- 
tion 0(£i,/f) = J c/ 2 k 1T 0(^!, k 1T ,fi) is the 6 — > limit of the initial condition 
0(£i,6, /i) of the Sudakov evolution. This approximation is reasonable be- 
cause of the strong suppression of / at large 6. Obviously, / is independent 
of the vector n. Now we make n approach t>2 (the light cone), and 0(£i,/z) 
coincides with the standard distribution function with the gauge invariance. 

The resummation for the distribution function of hadron 2 gives a simi- 
lar Sudakov factor. Hence, the additional smearing factor compared to the 
standard NLO QCD calculations is given by, 

s{ti,&,s,b) = mi P i,b)f(^,b) . (i6) 

Transformation of Eq. @ back to the kx space with Eqs. (14) and (|i"6|) 
inserted leads to 

da(s,p T ) 



dpi 



x#(£i, 6, s, p' T )S(£u&i s, k T ) 



(17) 
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with the variable p' T = |pr — k^|. We have set \i of (j) and H to the char- 
acteristic scale p' T of H. That is, the distribution functions <fi will evolve 
to p' T according to the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi equation 
[ |T6|j , which further incorporates the summation of the single logarithms \np' T . 



Equation (\T7\) is the factorization formula for direct photon production with 
the single k? approximation. 

To be related to the analysis in Q, we shall not work on Eq. dl7|) , 
but neglect the £ dependence of S, making the approximation S(s, kx) ~ 
5(0.5, 0.5, s, kr)- It is reasonable to assume £1, £2 ~ 0.5. Performing the 
integration over £1 and £2, Eq. (|ITD reduces to 



da(s, PT )_ r# kT MsM s{Sihrh (18) 



dpT J dp'q 
where the differential cross section 
da(s,p' T ) 



dp'rp 



J d^^M^MH^^sM . (19) 



is identified as the standard NLO QCD predictions, and S(s, kr) as the 
smearing function we shall investigate. Note that S for Figs. 1(a) and 1(b) are 
different, and thus Eq. (fHf ) in fact represents the sum over the two diagrams. 
While the smearing function employed in [Ij] is the same for Figs. 1(a) and 
1(b). 



III. NUMERICAL ANALYSIS 

Before proceeding with the numerical analysis of Eq. (|T8|) , we examine 
the smearing effect of the Sudakov factor S. It is known that the Gaussian 
smearing function exp(— k^/T 2 ) employed in |l| for direct photon production 
possesses the s-dependent width T, which is summarized as 

T < 1 GeV for ~ 30 GeV , 

r ~ 1 - 2 GeV for ^ ~ 30 - 100 GeV , 

r ~ 2 - 3 GeV for ~ 100 - 600 GeV , 

r ~ 3 - 4 GeV for ~ 600 - 1800 GeV . (20) 
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That is, T increases with s slightly. The predictions are very sensitive to the 
s dependence of T. Hence, it is nontrivial that the experimental data can be 
explained by our formalism. 

We demonstrate that the factor S(s, kx) exhibits the desired behavior 
shown in Eq. (|20|). For simplicity, we parametrize /(£^/s/2, b) as a Gaussian 
function exp(— T 2 b 2 /A) with T 2 = c(£i/i) r , which is the Fourier transforma- 
tion of exp(— fc^/r 2 ). The momentum fraction £ will be set to 1/2 in the 
numerical analysis below. The constants c and r are determined from the 
best fit of the parametrization to /, considering the variation of b between 
and 1/Aqcd = 5 GeV -1 , and of \fs between 20 and 1800 GeV. The re- 
sults are c = 0.12 and r = 0.60, if the parton is a quark, and c = 0.09 and 
r = 0.83, if the parton is a gluon. The decrease of f(Js/2, b) with b for dif- 
ferent y/s and its corresponding parametrization exp(— C\fsb 2 j^) from best 
fit are shown in Fig. 4. The values / = 1 at small b come from the truncation 
of / > 1 argued before. It is observed that the s dependence of the width T 
is roughly in agreement with Eq. fl20|). 

We then compute the differential cross section of direct photon production 
based on Eq. fll8"|) with the parametrized smearing function, 

5( S ,A; T )=exp[-4/(r 2 + r 2 )], (21) 

where the widths T^) = c(y / s/2) r are associated the hadron 1(2). We can 
certainly employ the exact Sudakov factor /, and convert it into the kx space 
numerically, when evaluating Eq. fll^l). This approach is not followed here, 
simply because we intend to make our analysis in full analogy with that in [p]]. 
In Fig. 5 we show the deviation (Data -Theory)/Theory of the NLO QCD 



predictions, obtained using the CTEQ4M parton distributions [[17]] , from the 
experimental data as a function of x t = 2pr/ y/s. Obviously, the deviation is 
huge, especially at low x t of each set of the data. In Fig. 6 the theoretical 
predictions come from Eq. (|i~8[ ) which includes the kr smearing in Eq. (|21~1) . 
It is clear that a significant improvement on the agreement between theory 
and experiments is achieved. 



IV. CONCLUSION 

In this paper we have identified the k? smearing, which is essential for 
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the explanation of the direct photon production data, as the Sudakov factor 
from the resummation of large radiative corrections to the parton distribution 
functions. The identification is confirmed by examining the s dependence of 
the Gaussian widths of the parametrized Sudakov factors, and the consis- 
tency of the improved predictions with the data. Compared to our analysis, 
the smearing employed in |1| is very naive. The effects should depend on the 
type of partons, and on the momentum fractions. A more accurate analysis 
even involves the recalculation of the hard scattering amplitudes with the 
partons off shell by k\x and which are associated with the hadron 1 
and 2, respectively. Our formalism including the resummation at low kx can 
be generalized to other QCD processes. All these subjects will be studied 
elsewhere. 

This work was supported by the National Science Council of R.O.C. under 
Grant Nos. NSC87-2112-M006-018 and NSC87-2112-M007-040. We also 
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Figure Captions 

Fig. 1. (a) and (b) Lowest-order diagrams of direct photon production. 
(c)-(h) 0(a s ) corrections to (a). 

Fig. 2. Factorization of direct photon production from hadron collisions. 

Fig. 3. (a) The derivative pfd^/dpf in the axial gauge, (b) General 
diagrams for the functions K and G. (c) The 0(a s ) function K. (d) The 
0(a s ) function G. 

Fig. 4. (a) Dependence of /(^/s/2, b) on b (dashed line) and its correspond- 
ing parametrization (solid line) for (1) y/s = 30 GeV, (2) y/s = 600 GeV, 
and (3) y/s = 1800 GeV associated with (a) the quark distribution function 
and (b) the gluon distribution function. 

Fig. 5. Compilation of direct photon experiments compared to the NLO 
QCD predictions using the CTEQ4M parton distributions 

Fig. 6. Compilation of direct photon experiments compared to the kr- 
resummed predictions using the CTEQ4M parton distributions 
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